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ABSTRACT 

The hierarchical clustering model of the cosmic mass field is examined 
in the context of intermittency. We show that the mass field satisfying the 
correlation hierarchy £ n ~ Qnfe)™ -1 is intermittent if k < d, where d is the 
dimension of the field, and k is the power-law index of the non-linear power 
spectrum in the discrete wavelet transform (DWT) representation. We also find 
that a field with singular clustering can be described by hierarchical clustering 
models with scale-dependent coefficients Q n and that this scale-dependence 
is completely determined by the intermittent exponent and k. Moreover, the 
singular exponents of a field can be calculated by the asymptotic behavior of 
Q n when n is large. Applying this result to the transmitted flux of HS1700 
Lya forests, we find that the underlying mass field of the Lya forests is 
significantly intermittent. On physical scales less than about 2.0 h _1 Mpc, the 
observed intermittent behavior is qualitatively different from the prediction of 
the hierarchical clustering with constant Q n . The observations, however, do 
show the existence of an asymptotic value for the singular exponents. Therefore, 
the mass field can be described by the hierarchical clustering model with 
scale-dependent Q n . The singular exponent indicates that the cosmic mass field 
at redshift ~ 2 is weakly singular at least on physical scales as small as 10 h _1 
kpc. 

Subject headings: cosmology: theory - large-scale structure of universe 
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1. Introduction 

The cosmic mass density field seems to posses two basic features which at first glance 
appear contradictory. The first is the statistical homogeneity and isotropy of the field. The 
second is that the mass field consists of isolated density peaks like galaxies and clusters of 
galaxies, in which the mass densities are much higher than average. Moreover, coherent 
structures, like filaments, sheets and even cellular or network structures with characteristic 
lengths have been detected in galaxy distributions. 

These two features can be reconciled by assuming the formation of a highly intermittent 
field from an initially gaussian random field. In this scenario, the growth of structure in 
the mass field is initially linear and randomness is the mechanism by which structure is 
formed. At later stages, non-linearities enter the dynamics and prevent infinite growth. 
As a result, the structures that thereby arise in the random mass field have an interesting 
character: they are strong enhancements (peaks) of the density field scattered in a space 
with a low density background. This feature generally is called intermittency and was 
originally introduced for describing the temperature and velocity distributions in turbulence 
(Batchelor & Townsend 1949; for review, see Frisch, 1995.) An early attempt of studying 
the formation of the large scale coherent structures from a random field is the theory of 
pancakes (Zeldovich 1970). This work is the predecessor of the intermittency approach 
to the formation of coherent structures (Shandarin & Zeldovich 1984; for a review, see 
Shandarin & Zeldovich 1989). 

The intermittency picture is supported by modelling the cosmic matter distribution 
in non-linear regime with lognormal(LN) random field, which is a simplest model of 
intermittency behavior in random fields (e.g. Zeldovich, Ruzmaikin, & Sokoloff 1990). The 
LN model has been found to reporoduce the nonlinear matter distribution evolved from 
Gaussian initial conditions and gives rise to spatial patterns characterized statistically by 
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a kind of intermittency ('heterotopicity') (Cole and Jones, 1991; Jones, Coles & Martinez, 
1992). The LN model of the baryonic matter distribution (Bi & Davidsen, 1997) has also 
been successfully applied to explain the non-Gaussian features of the underlying mass field 
traced by Lya forests (Feng & Fang 2000). 

The intermittency picture is also supported by the singular density profile of massive 
halos (like p{r) oc r~ 7 with exponent 7 ~ 1) given by N-body simulations (e.g. Navarro, 
Frenk & White, 1997, Jing & Suto 2000.) This profile indicates that the events of large 
differences between the local densities separated by small distances r are more frequent than 
with a gaussian field. As such the probability distribution function(PDF) of the density 
difference, A r (x) = p(x + r) — p(x), on small scales r is long tailed. 

The PDF of A r (x) can only be effectively detected by a scale-space decomposition 
because A r (x) consists of both scale r and position x. A scale-space decomposition analysis 
of the mass density field traced by the Lya forests has indeed obtained more direct evidence 
of the intermittency (Jamkhedkar, Zhan, & Fang 2000.) They found that the PDF of 
A r (x) is gaussian on scale r larger than about 2 h" 1 Mpc, but becomes longed tailed as 
r decreases. Hence, the random mass field traced by the Lya transmitted flux shows an 
excess of large fluctuation on small scales in comparison to a Gaussian distribution. This is 
typical of an intermittent field (Shraiman & Siggia 2000.) 

It seems that intermittency is a common feature of the non-linearly evolved cosmic 
mass field. At the same time, hierarchical clustering is believed to be a common feature 
of the non-linearly evolved cosmic mass field, and has been widely applied to construct 
semi-analytic models of gravitational clustering in the universe (e.g. White 1979). Although 
hierarchical clustered fields show indeed some kind of intermittency, it is unclear whether 
the correlation hierarchy is consistent with the observed intermittency of the cosmic mass 
field. For instance, in the LN model, the high order correlations do not satisfy a hierarchical 
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relation, but obey a well-known Kirkwood scaling relation (Kirkwood, 1935; Peebles, 1980). 
Therefore, a LN random field is intermittent, but does not satisfy the hierarchical clustering. 

The purpose of this paper is to investigate the relationship between intermittency and 
hierarchical clustering. We will not limit ourselves to special models of intermittency, but 
will try to calculate the intermittent feature of a hierarchical clustered field in general. 
We first introduce the intermittent exponent which gives a detailed classification of the 
non-linear clustering. Specifically it gives a complete and uniform description of the 
nonlinearity of the cosmic mass field from initially gaussian perturbations, to weak, 
and eventually strong non-linear clustering. We then study the intermittent exponent 
predicted by the hierarchical clustering. As expected, the intermittent behavior given by 
the hierarchical clustering model is found to be in good agreement with observed results on 
large scales. However, for the highly non-linearly evolved field, we find that the simplest 
correlation hierarchy, i.e., one with constant coefficients Q n , is no longer adequate to 
describe the observed features in the density field. In this case, the model of hierarchical 
clustering with scale dependent Q n may still work. Moreover, we find the relationship 
between the intermittent and singular behaviors of the cosmic mass field and the scale 
dependence of the Q n . 

The paper will be organized as follows. §2 introduces the basic statistical measures of 
the intermittency of random mass fields - the structure function and intermittent exponent. 
§3 presents the predicted intermittency of a hierarchically clustered field. The singular 
behavior of the hierarchical clustering with scale-dependent Q n will be discussed in §4. The 
comparison of the intermittency and singular behaviors of a hierarchically clustered field 
with samples of QSO's Ly-a forests is given in §5. Finally, the conclusions and discussions 
are in 86. 
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2. Intermittency 

The chief characteristic of an intermittent field is that the interesting stuff happens 
in a local area. As a consequence, the description of intermittent fields by traditional 
methods such as the Fourier power spectrum are ineffective because these methods are 
essentially non-local. For instance, the amplitudes of the Fourier coefficients lose the spatial 
information. Two fields which have a similar gaussian PDF of the amplitudes of the Fourier 
coefficient may have very different intermittent features. Though the phases of the Fourier 
coefficients do contain the spatial information, they are not computationally convenient. 

The density profiles described in §1 look to be a good way to the measure mass density 
peaks. The profiles are, however, only a measure for the individual massive halos and not 
a statistical measure of the entire random field consisting of these peaks. For example, 
imagine a homogeneous gaussian field in a cubic box with 128 3 grids, which corresponds to 
f 28 3 realizations of a gaussian random variable. Rare peaks as high as 4<r will occur at some 
points. Obviously, the individual density profiles are useless in measuring the excess of the 
large density fluctuations with respect to a gaussian distribution. Moreover, the exponent 
7 introduced with the individual density profiles p(r) oc r -7 cannot be used to describe 
a statistically homogeneous and isotropic field, which requires that proper measures be 
statistically invariant quantities under translational and rotational transformations. 

The proper measure of an intermittent field should satisfy the following conditions 
simultaneously: f . its second order statistics should describe the power spectrum of the 
field, and 2. the correlations of phases should describe coherent structures. We show in this 
section that the structure functions in the discrete wavelet transform (DWT) representation 
meet these conditions. 
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2.1. The local density difference 

The basic quantity for describing the intermittency of a random density field p(x) is 
the difference between the densities at the position x with separation r, i.e., 

A r (x) = p(x + r) - p(x) = l[*(x + r) - 5(x)], (1) 

where <5(x) = [p(x) — p]/p is the density contrast. For simplicity we assume that the density 
is normalized, p — 1. The density difference is essential to describe the most important 
intermittent characteristics of a mass field. These features include: 

f . The power spectrum and the local power spectrum. 

The ensemble average of the second moment of A r (x) is 

S 2 (r) = (|A r (x)| 2 >. (2) 

If the field is homogeneous, S 2 (r) is independent of x and depends only on r. S 2 (r) is the 
mean power of the density fluctuations at wavenumber k ~ 2n/r, and therefore, S 2 (r) is 
the power spectrum of the field (see §2.3). For a given realization of the random field, 
|A r (x)| 2 is the local power spectrum at x. For an intermittent field the spatial distribution 
of |A r (x)| 2 is highly irregular or spiky. The power of the density fluctuations on small r is 
highly localized in space with very low powers between these spikes. This is probably the 
easiest way of identifying intermittency (Jamkhedkar, Zhan, & Fang 2000.) 

2. The long-tailed PDF of A r (x). 

One of the defining characteristics of an intermittent field is the long-tailed PDF of 
A r (x) on small scales r. This trait can effectively be measured by the higher order moments 
of A r (x), 

S 2 » = (|A r (x)| 2 ™> (3) 
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where n is a positive integer. The S 2n (r) are called structure functions. When the "fair 
sample hypothesis" is applicable (Peebles 1980), S 2n (r) can be calculated by the spatial 
average, 

S 2n {r) = ^J |A r (x)| 2 "dx, (4) 

where V is the spatial normalization. When n is large, S 2n (r) is dominated by the long-tail 
events. Eqs. (2) and (4) show that the structure function S 2n unifies the analysis of 
intermittency (n > 1) with the power spectrum (n = 1). 

3. Density peaks and scale-scale correlations. 

The density profile of mass halos p(|x — xo|) oc |x — xo|~ 7 at xo means that events 
with large |A r (x)| on different scales r occur at the same place x = x . In other words, the 
density peaks lead to the correlation between |A r (x)| and |A r /(x)| with r ^ r'. As such, 
one can distinguish the field containing only peaks caused by gaussian fluctuations from an 

intermittent field by the scale-scale correlations defined as (Pando et al. 1998) 

, (A r (x)»A r ,(x)"') 
r > r ' (A r (x)«)(A r ,(x)"')' 1 ) 

where n and n! are even integers. C™y = 1 for gaussian fields, while C™y > 1 for 
intermittent fields. When n and n' are large, C™'™ is dominated by the high density 
peaks. The excess of high peaks (as compared to gaussian peaks) can be measured by the 
scale-scale correlations. 

4. Singular behavior and Holder exponent a 

In the context of random fields, the singular behavior is characterized by the Holder 
exponent a defined as (Adler 1981) 

A r (x) = |5(x + r) -5(x)| < (r/L) a as r^O, (6) 

where the constant L can be taken as the sample size. The exponent a measures the 
smoothness of the field: for larger a, the field is smoother on smaller scales and vice versa. 
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If a is negative, the field becomes singular. In this case the Holder exponent is dominated 
by the events with extremely large density difference |A r (x)| on extremely small scales 
(r — > 0), and therefore, the exponent a should be dependent on the index 7 of the singular 
density profile if the object with the singular density profile is typical of the mass field. 
For a homogeneous and isotropic field, the PDF of the local density (contrast) difference 
[5(x + r) — <5(x)] for a given scale r has to be independent of position x. Therefore, the 
exponent a defined by eq.(6) is statistically invariant under translation and rotational. 

Since the resolution of real and simulation samples is always finite, one cannot 
practically measure the individual singularity of |x + r| -7 at position x with r ~ 0. 
However, a defined by eq.(6) can be calculated as the asymptotic behavior of the statistics 
of local density differences when 1/r — > 00. Hence, one can measure the singular exponent 
a as a statistical property of the random field, rather than using the density profiles of the 
individual massive halos. 

We should make an important distinction between statistics based on A r (x) and 
statistics that use the density p(x) or density contrast 5(x). The former is the difference 
in either the density p(x) or the density contrast <5(x) at different positions in space. It 
is possible that a large difference A r (x) on small scales occurs in regions of low density, 
and also possible that a small difference occurs in regions of high density. To study 
intermittency, data which samples the entire cosmic mass field is necessary. 

Finally, an advantage of using the local density difference is in reducing the 
contamination of the random velocity field. It is well known that the random velocity 
field will smooth the mass field in redshift space and repress the power of the density 
perturbations on scales below that characterized by the velocity dispersion. Because the 
difference A r (x) is localized in space and scale, the velocity fluctuations on scales larger 
than r will lead to offset of the event A r (x) in redshift space, but keep the magnitude of 
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A r (x) unchanged. As a consequence, the PDF of A r (x) will not be affected by velocity 
fluctuations on scales larger than r, and the spiky features of the local power spectrum of 
the Lya forests are still significant on scales of a few tens h -1 kpc (Jamkhedkar, Zhan, & 
Fang 2000.) 



2.2. The Intermittent exponent 

In the previous section we showed that the intermittency is described by the 
scale-dependence of the higher order moments of A r (x). In this section, we introduce the 
intermittent exponent, which is a more effective tool of describing this scale-dependence. 

For a homogeneous and isotropic gaussian field, the structure function eq.(3) can be 
calculated by 



S 2n (r) = / P g (A r )[A r } 2n dA r (7) 

J — oo 

= (2n-l)!![^(r)]", 

where P g (A r ) is the gaussian PDF of A r (x), and a 2 (r) = (|A r .| 2 ) is the variance of P g (A r ). 

The deviation of the field from a gaussian distribution on scale r can be measured by 
the intermittent exponent defined asQ 

( is equal to zero if the mass density field is gaussian on all scales. From eq.(8), we have 



S 2n (r ) \r 



4 In turbulence, S 2n (r) is used to define the so-called anomalous scaling describing the 
intermittent behavior (Shraiman & Siggia 2000.) 
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where again, ( is equal to zero if the mass density field is gaussian on all scales. 
A field is self-similar if the difference A r (x) as a random variable satisfies 

A r (x) = A A A Ar (x), (10) 

where h is constant. In this case, S 2n (r) is independent of r. Once again, the exponent ( 
defined by eq.(9) is zero. A field is said to be intermittent if the exponent ( is non-zero. 
Hence, an intermittent field is neither gaussian nor selfsimilar (Frisch, 1995.) 

If the field contains more "abnormal" events, i.e., large A r (x) on small scales, S 2n (r) 
will be dominated by these events, especially when n is large. In this case, we have 
S 2n (r) > [S 2 (r)] n on small r, and therefore, the exponent ( is less than zero. Similar to the 
exponent a [eq.(6)], the intermittent exponent ( measures the smoothness of the field: for 
larger (, the field is smoother on smaller scales, and vice versa. If ( is negative on small 
scales, the field is probably singular. 

The n- and r-dependence of ( provide a detailed description of the non-linear features 
of the mass field. ( gives an unified criterion for classifying the non-linear features of a 
random field, from gaussian, to self-similar, to mono- and multi-fractal, and to singular. 

2.3. The Intermittent exponent in the DWT representation 

The density contrast difference A r (x) [eq.(l)] contains two variables: the position x 
and the scale r, and therefore, the information of A r (x) can best be extracted by a proper 
space-scale decomposition. 

We will use the discrete wavelet transform (DWT) decomposition [For details on 
the DWT refer to Mallat (1989a,b); Meyer (1992); Daubechies, (1992), and for physical 
applications, refer to Fang & Thews (1998)]. To simplify the notation, we consider only 
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the DWT decomposition for 1-D density field p{x) or S(x) on a spatial range L. It is 
straightforward to generalize the result to 2- and 3-D fields. 

We first divide the spatial range L into 2 j segments labelled by I — 0, 1..2 j — 1. Each 
of the segments has length L/2 j . The density contrast difference A r (x) = S(x + r) — 
defined in eq.(l) can be replaced by 

"(J+l/2)L/2>' r (l+l)L/2i 



e 3,l = \ T 



IL/23 



S(x)dx — / 
^0 



(Z+1/2)L/2J 



5{x)dx 



(11) 



e J: ; measures the difference between the mean density contrasts in the localized segments 
lL/2 j < x < (I + l/2)L/2 j and (I + l/2)L/2 j < x < (I + l)L/2 j . Therefore, it can be 
identified as A r (x) with x = lL/2 j and r = L/2 j . 

is the Haar wavelet function coefficient (WFC), given by the projection of S(x) onto 
the Haar wavelet basis ^i(x) as 



= J 5(x)^f :l (x)dx, 



(12) 



where 



1 from Ll2~i to Ll2~i + L2^- 1 

-1 from L12~^ + L2^- x to L(l + l)2"-» 

otherwise. 



(13) 



The factor y2i/L insures the orthonormality of the wavelet with respect to both indices j 
and 

For other wavelets vpj : i(x), the WFC, e,y, also measures the difference of the density 
contrast on separation L/2 J at a position lL/2 j . We will use the Daubechies 4 wavelet 
(Daubechies, 1992) in our numerical calculations below because it is better behaved in scale 
space than the Haar wavelet (Fang & Thews 1998). The set of wavelets ipj,i(x) (j = 0, 1, .. 
and / = 0, 1, ..2 j — 1) is complete and orthogonal. In the DWT representation, the field 5(x) 
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can be expressed as 

J 2>'-l 

= E E Zi'rti'A*)- ( 14 ) 

j'=o ;=o 

J is determined by the spatial resolution Ax of the sample, i.e., J ~ log 2 (L/Aa;). The 
WFC e,',/ is given by eq.(12), replacing the Haar wavelet i])fi{x) by ipj t i(x). Since wavelets 
are admissible, the WFC can be calculated with either the density or the density contrast 

= J p{x)^j,i{x)dx = j=/ S(x)ip jt i(x)dx. (15) 
In the DWT representation, eq.(2) is 

^ = ^eW- (16) 

This is the power spectrum in the DWT representation and equal to a band-averaged 
Fourier power spectrum given by (Fang & Feng 2000) 

1 oo 

S ) = T 3 E \4>(n/n 2 P(k), (17) 

n=— oo 

where -P(fc) is the Fourier power spectrum with the wavenumber k = 1-nnjL. ip(n) is the 
Fourier transform of the basic wavelet ip(i]). 

The structure function (3) in the DWT representation is 

sf HI^ 2n ) = ^EI^I 2n (18) 

1 1=0 

and eq. (8) is 

ST = j^r - (1 9 ) 

It has been shown that ( given by eq.(19) is the same as (8) (Jaffred 1994.) S| n or ( is the 
basic measure of the intermittency of the cosmic mass fields. 
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3. Intermittency within hierarchical clustering 

Hierarchical clustering provides a model describing the non-linear clustering of the 
cosmic mass field. It assumes that the correlation functions of the mass density can be 
described by the linked-pair approximation, i.e., the n-th irreducible correlation function £ n 
is given by the two-point correlation function £ 2 as £ n = Qn^' 1 ? where Q n is the hierarchical 
coefficient (White 1979). As the high order correlation functions provide a direct way to 
investigate the nonlinearity in a random field and test the scaling hierarchy, it follows that 
one can study the hierarchical clustering predictions pertaining to intermittency and the 
singular behavior of the mass field. In this section we develop the higher order correlation 
functions in the DWT representation. 



3.1. Linked-pair approximation with the DWT modes 

For the n-th order correlation function, the correlation hierarchy, or linked pair 
approximation, is 

Uxi,..x n ) = E«E S^W, (20) 

™ (ab) 

where coefficients are constant, the sum w is over the types of the tree graphs with n 
vertices, x a b = |x a — x b |, and the sum (ab) is over relabelings within w (Fry 1984). With 
this relation, the higher order behavior of the clustering is completely determined by the 
two-point correlation £2(^06 ) and the constants Q™- Hence, the intermittent behavior of 
a hierarchical clustered field can be determined by the coefficients and the non-linear 
power spectrum. 

For instance, in the case of n — 4, eq.(20) yields 

(5( Xl )5(x 2 )5(x 3 )5(x 4 )) = (21) 
Q2[(5(x 1 )5(x 2 ))(5(x 2 )5(x 3 ))(5(x 3 )5(x4)) + eye. 11 terms] 
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+g^[<5(x 1 )5(x 2 ))(5(x 1 )5(x3))(5(x 1 )5(x 4 )) + eye. 3 terms]. 

where Q% is for snake diagrams, and Q\ is for stars. 

We express eq.(21) in the DWT basis, ^j4(xi)^j i i(x 2 )V',j,i(x3)V ; j,i(x4), where ^j,i(x) is 
the 3-D wavelet (x = x 1 , x 2 , x 3 ). ip$,\(x) is given by a direct product of the 1-D wavelets, 
i.e., Vj,i(x) = $ jl j 1 (x 1 )rl; jath (a?)i(> j3ih (x 3 ). 

Equation (21) yields 

= QtiY, E S ( g j,i g j',i')( g j",i" g j"',i'")(ej"",i""ej,i> (22) 

j/ j/ j// j// j/// yrr yrrt yttt 
J ^j,l(x 2 )^ j / iF (x 2 )^ j » i i//(x 2 )rfx2 y V'j,l( x 3)V'j"',l"'( X 3)V'j"",l""( x 3)rfX3 

+cyc. 11 terms] 

E (ejM'ej,i)<ej" ) i»ej 1 i)(ej/« )1 /«ej i i) 

j/ y y f \" y ,f \ fff 

J V'j,i(xi)V'j',i'(xi)V'j'',i"( x i)V'j'",i'"( x i) rfx i +cyc- 3 terms] 
where the j and 1 terms are the 3-D scales and positions respectively. 

The DWT is very efficient in compressing data, i.e., the off j-diagonal elements of the 
covariance generally are much less than the j-diagonal elements. This property is 

also consistent with the quasi-locality of mass clustering (Pando, Feng & Fang 2001.) Thus, 
the r.h.s. of eq.(22) is dominated by the terms (e?i) 3 and we have 

{ll)^\QlA- i + Q\B^)\ (23) 

where the factors Aj and Bj are given by 



4 



/Vf,i(x)dx] 2 (24) 
/ ^i( x ^ x - (25) 
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Because wavelets are constructed by dilating and translating a basic wavelet ip(rj), A- s 
and B- } can be more simply written as 

l2 



A. - of+f+f J_ 
J L 3 



tp {r])drj 



and 



(26) 



(27) 



Thus, eq.(23) becomes finally 



(^)^Q 4 2^- 3 (^) 3 , 

where the coefficient Q4 is independent of j and given by 

f 



(28) 



Qa 



L 3 



(29) 



Similarly, for 2n-th order, we have 



(eff) - Q 2 n2 {n_1)(il+i2+i3) (e? 1 ) 2n " 1 , 



(30) 



where again Q 2n is independent of j. Eq.(30) is the counterpart of correlation hierarchy 
£n = Qni2~ l m the DWT representation. 

For a 1-D field, such as Lya forests, we use a projection of a 3-D distribution <5(x) onto 
1-D as 

/oo 
5(x)ip j j(x 1 )(f)j^ m (x 2 )(f)j< tn (x 3 )dx 1 dx 2 dx 3 , (31) 
-00 

where x 1 is the redshift , x 2 and are the positions in the sky, and 4>j,i{x) is the scaling 
function of the DWT analysis. The scaling function (f>j,i(x) plays the role of a window 
function on scale j at position I (Fang and Thews 1998). With eq.(31), the 1-D field 
linked-pair approximation (20) yields 



(6j7)^g 2n 2^"- 1 )(^) 2 - 1 , 



(32) 
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where Q 2n is for a 1-D field and is generally not equal to Q 2n for a 3-D field because the 
difference in the projection (31) for the 1-D case. We use the same notation Q 2n for both 
1-D and 3-D and from the context it will be clear whether we are using the 1-D or 3-D 
quantity. 



3.2. The intermittency exponent within hierarchical clustering 

Applying eq.(30) for the diagonal case, i.e. j 1 = j 2 = j 3 = j, we have 

(eS)^g 2 n2 3j( "- 1) (6j il ) 2 ^ 1 . (33) 

One can rewrite eqs.(32) and (33) as 

(6fP^g 2 „2* ( "- 1 )(^ 1 ) 2 "- 1 , (34) 

where d is the number of dimensions. 

Substituting eq.(34) into eq.(19), we have 

S? n = Q 2n 2 d ^ n - 1 \iD n - 1 . (35) 

Therefore, the intermittent exponent of a hierarchical clustered field is 

Ch = -(n-l) 



d + - ln 2 Pj 

J 



--ln 2 Q 2n , (36) 



J 

the subscript h is used to emphasize that (h is the predication in the hierarchical clustering 
model. For small scales (large j), the last term of eq.(36) is negligible if the Q 2 nS are 
constant. Hence, for a hierarchically clustered field, the intermittency exponent (h is 
completely determined by the power spectrum Pj = (e 2 ^) in the non-linear regime. If the 
power spectrum is a power law, Pj oc 2~ JK , we have 

C h ~-(n-l)(d-«). (37) 
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Hence, the field is intermittent if k < d. 

The intermittent exponent eq.(37) is independent of j. This is the simplest type of 
intermittency - monofractal with fractal dimension k. That is, the space filling of the 
density perturbations will be less by a factor (l/2)( d-/t ) from scale j to j + 1, and in effect, 
the hierarchical clustering model with constant Q n is equivalent to assuming that the 
random field 5(x) is monofractal. This conclusion is consistent with the phenomeno logical 
model of the hierarchical relations developed by Soneira & Peebles (1977). The model is 
essentially based on a dimension k fractal distribution.^ 



4. Singular behavior of hierarchical clustering 

4.1. Singular exponent 

We now consider the singular exponent a of eq.(6). If the mass field contains singular 
structures like p oc r~ 7 with 7 > at position x (or 1 ), the modulus of the diagonal WFCs 
(ji = 32 = h = j) satisfy 

<A2^~ d ' 2 \ (38) 



where the factor d/2 is from the normalization y2 J /L in eq.(13). Eq.(38) shows that the 
singular exponent defined in eq.(6) is a = 7 — d/2. Thus, one can measure 7 asymptotically 
as follows. 

Because for large n, S 2n is dominated by events with large e^, eq.(19) gives 

7 x C, H ; ln 2 Pj, when n is large. (39) 

2 2 tx — I 



5 This model is also similar to the so-called /3-model of the intermittency of turbulence, 
for which the exponent ( has similar terms as eq.(42) (Frisch 1995.) 
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Thus, a field is singular if the following condition holds 

- — — ( + — ln 2 Pj > 0, when n is large. 



(40) 



For a hierarchical clustered field, condition (40) is 



d — k > when n is large. 



(41) 



Therefore, a hierarchical clustered field cannot be singular if the power index k > d. Under 
condition (41), eq.(37) leads to ( < 0. This result is consistent with the statement that the 
field is singular when ( < on small scales (§2.2). 



In order to have the hierarchical clustering scenario compatible with observed data 
and N-body simulation samples, it is often assumed that the Q n are scale dependent. In 
this section, we will show that the scale dependence of the Q n is given by the intermittent 
exponents. The validity of the hierarchical clustering model with scale dependent Q n can 
be tested by examining the relationship between the Q n scale dependence and the singular 
exponents. 

For the hierarchical clustering models with scale-dependent Q2 n , eq.(34) still holds, as 
all the quantities in eq.(34) are on same scale j. Thus, the scale-dependence of Qm is given 



4.2. Hierarchical clustering with scale dependent coefficients 




(42) 



The scale dependence of Q 2n can be measured by a power law index defined by 



Q 2 n OC 2^ 2 ". 



(43) 



Using eqs.(19) and (42), eq.(43) yields 



fan = ~d(n - 1) - C 



n- 1 



ln 2 P r 



(44) 



J 
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We see that the scale-dependence of Q 2n is completely determined by the intermittent 
exponent ( and the non-linear power spectrum. 

Using eq.(39), the singular exponent 7 can be expressed by (5 2n as 

1 1 

7xd+ — f3 2n + - ln 2 Pj, when n is large. (45) 
2n j 

Eq. (45) shows that the singular index 7 is determined by the scale-dependence of the 
coefficients Q 2n and the non-linear power spectrum Pj. 

The important point of eq.(45) is that the l.h.s. is given by the singular profiles such 
as p oc r~ 7 , independent of n, while the r.h.s. consists of n-dependent quantities like f3 2n . 
Therefore, a test of the scale-dependent Q n model is to check whether 7 given by eq.(45) is 
n-independent when n is large. 

It is known that a scale-dependent intermittent exponent indicates that the mass field is 
no longer monofractal, but is instead, multifractal (Farge at al. 1996.) Galaxy distributions 
have been shown to be multi-fractal in nature (Jones et.al, 1988) and to explain this 
distribution, a phenomenological model of non-linear fragmentation with multi-scaling 
parameters has been proposed (Jones, Coles & Martinez, 1992). Our present study has 
revealed the intrinsic relationships between the scale-dependence of Q 2n , multifractal 
nature, and singular behavior of the mass field. 

This completes our derivation of intermittency within hierarchical clustering. With 
eq.(18) we pick up the DWT power spectrum, eq.(37) gives the intermittency exponent 
predicted by the hierarchical picture, eq.(42) provides a way to check for the scale 
dependence of the correlation coefficients, and eq.(45) gives a measure of singular behavior 
of the field. The versatility of this approach is clear. 



5. Intermittent exponents of Lya forests 



Since the intermittent and singular features of the cosmic mass field are measured by 
the scale-dependence of the structure functions, samples covering a large range of scales are 
best suited for this analysis. As has been emphasized, the intermittent and singular features 
are traits of the mass field and not of the individual massive halos. As such, the transmitted 
flux of QSO's Lya forests given by high resolution absorption spectrum are ideal data sets 
with which to work. The Lya transmitted flux is due to absorption by gases in cool and 
low density regions. The pressure gradients are generally less than the gravitational forces. 
The distribution of cool baryonic diffuse matter is almost point-by-point proportional to 
the underlying dark matter density (Bi, Ge & Fang 1995). The statistical features of the 
underlying mass field of the Lya forests can be detected by the transmitted flux. 



5.1. Samples 

We use the Lya transmitted flux of QSO HS1700+64 (z = 2.72) for our analysis. This 
sample has been employed to study the baryonic matter density (Bi & Davidsen 1997), 
the Fourier and DWT power spectra and non-gaussian features (Feng & Fang 2000). The 
recovered power spectrum has been found to be consistent with the CDM cosmogony on 
scales larger than about 0.1 hr x Mpc. The data ranges from 3727.012A to 5523. 554A, for a 
total of 55882 pixels. On average, a pixel is about 0.028A, equivalent to a physical distance 
of ~ 5 h _1 kpc at z ~ 2 for an Einstein-de Sitter universe. These scales weakly depend on 
cosmological density parameter Q. One can safely ignore the effect of the non-linear relation 
between redshift and distance for our present analysis. In this paper, we use the data from 
A = 3815. 6A to 4434. 3A, which corresponds to z = 2.14 ~ 2.65. The comoving spatial size 
of our data is 189.84 h -1 Mpc in the CDM model. This means for scale j, the spatial size is 
189.84/2^' h" 1 Mpc. 
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With Lya data there is always the danger of contamination from the presence of metal 
lines. We use three ways to estimate the error this contamination causes. First, we block 
out the significant metal line regions. Since the WFC ijj are localized, the metal line 
regions are easily separated from the rest. Any effect from these regions can be removed 
by not counting the DWT modes in the blocked regions. The second way is to fill those 
regions with random data which has the same mean power as the rest of the original data 
and to smooth the data over the boundaries. Lastly we discard the metal line chunks and 
smoothly connect the good chunks of data. We find that the results not sensitive to the 
method of removing metal lines. 

Another source of contamination is noise. To estimate the effect of noise we smooth 
the QSO's spectrum by filtering out all extremely sharp spikes in the local power spectra 
on finest scales. These spikes are caused by relatively strong fluctuations between two 
neighboring pixels. Since such events are on the smallest scales only, the analysis on larger 
scales does not depend on whether we smooth the sample or not. 

The DWT power spectrum of the Lya transmitted flux of QSO HS 1700+64 is shown 
in Fig. 1. It is clear that k is larger than 1 on all small scales j > 9 or physical scales < 
100 h" 1 kpc. 



5.2. Observed intermittent exponent 

We calculate the intermittent exponent with eq.(19), 

C=— ln 2 Sf. (46) 
J 

The result is plotted in Fig. 2 and shows that on scales j > 9, the intermittent exponents 
are always significantly less than zero. The underlying mass field of the QSO forests is 
neither non-gaussian nor non-selfsimilar, but essentially intermittent on physical scales less 
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than 100 h" 1 kpc. 

We next test the predicted relation (36) between the intermittent exponent and power 
spectrum. The results of (h are plotted in Fig. 3. Because the power law index k of the 
DWT power spectrum is generally larger than 1, ( h is always positive. This is qualitatively 
different from ( on small scales. Therefore, only on large scales j < 5 or physical scales > 2 
h _1 Mpc, does the correlation hierarchy with constant Q n match observations. 



5.3. Multifractal scaling and singular exponent 

To check for hierarchical clustering with scale-dependent coefficients, we calculate Q 2n 
with eq.(42). The results are shown in Fig. 4, which shows that only for lower n, can 
the coefficient approximately be considered as constant. At all higher n, Qi n significantly 
increases with decreasing scale after j > 9, which implies the breaking of a single scaling 
hierarchy in the baryonic matter distribution. 

Fig. 5 plots the singular exponent 7 calculated by eq.(45). It indeed shows that 
for each j, 7 is asymptotically independent of n, when n is large. That is, models with 
scale-dependent Q n are reasonable and the asymptotic values of 7 are probably given by 
the exponent of the singular mass density profile. 

Fig. 5 also show that the singular exponents 7 become larger with increasing scales. 
For instance, 7 ~ 0.45 at j = 9, while 7 ~ 0.18 at j = 13. This is understandable if we 
consider the so-called "universal" density profile of massive halos (Navarro, Frenk & White 
1997) 

~ f I ^2 ( 47 ) 
r(r + a) z 

where a is the core radius. From this profile it is clear that 7 on small scale r < a will be 
less than for r > a. Comparing with (47), the 7 given by the Lya transmitted flux are small 
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and implies that the singular behavior of mass field at redshift z ~ 2 is much weaker than 
the present. 

Although Fig. 5 shows that 7 > on scale j = 13 or ~ 7 h _1 kpc, the field is not 
necessarily singular. For instance, if the profile (47) is modified as 

P(r) ~ (r + 6)(r + a)' (48) 
where physical scale b ~ 5 h _1 kpc, 7 will no longer be positive on scales j > 14. 



6. Conclusion and discussions 

The structure functions provide a complete and unified way of studying intermittent 
fields. With this method, we have shown that intermittency and hierarchical clustering in 
the non-linear mass field are related. We find that if the field is singular, the intermittent 
behavior on higher order n can be completely described by the singular exponents and 
the non-linear power spectrum. In this case, the scale-dependencies of the hierarchical 
clustering coefficients Q 2n at large n are also completely determined by these parameters. 

The transmitted flux of HS1700 does show the existence of the asymptotic value of the 
singular exponents. The cosmic mass field at redshift ~ 2 is probably weakly singular at 
least on physical scales as small as 10 h _1 kpc. A singular field can be neither self-similar 
nor mono-fractal. It is well known that under many different conditions the field given by 
gravitational clustering is self-similar. Therefore, the presence of singular clustering implies 
a broken self-similarity and/or mo no- fractal. 

From a dynamical point of view, there is a basis for understanding the relation between 
the intermittency and hierarchical clustering. The basic dynamical processes of hierarchical 
clustering are the merging of two halos or the accretion of small halos into bigger ones. 
Regardless of specific processes, either the infall of small halos into massive ones or the 
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merging of massive halos, these processes are stochastic. Thus, the hierarchical clustering is 
described by equations containing the addition and multiplication of stochastic variables. 
But these stochastic processes are typical of the origin of intermittency (e.g. Nakao 1998.). 
Even when the initial PDF of the random variables is gaussian, the dynamics involving 
additive and multiplicative stochastic variables will evolve to a long tailed PDF. For 
instance, the growth of the mass M of a halo is sometimes phenomenologically described 
by the rate equation of halo mass growth (e.g. Cole & Lacey 1996, Salvador-Sole, Solanes 
& Manrique 1998.) The rate equation is a phenomenological model of the formation of 
intermittency via merging. 

The so-called merging history tree theory has widely been employed to describe the 
mean merging history of individual halos and then to fit with the statistical properties of 
galaxies. Many merging tree models have been proposed. In the context of intermittency, 
different merging trees have different rules for merging from smaller to larger scales, and 
therefore, lead to different scale-scale correlations (§2.1). Thus, discrimination among 
merging tree models is possible via a intermittent analysis. 

We thank Dr. Wolung Lee for his support to this work. We also thank Dr. D. Tytler 
for kindly providing the data of the Keck spectrum HS 1700+64. LLF acknowledges support 
from the National Science Foundation of China (NSFC) and National Key Basic Research 
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Fig. 1.— The DWT power spectrum of the Lya transmitted flux of QSO HS1700+64. The 
top scale is the comoving scale, i.e., 189. 84/2-? h" 1 Mpc. 
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Fig. 2. — The intermittent exponent ( vs. j of the Lya transmitted flux of QSO HS1700+64. 
The top scale is the comoving scale, i.e., 189.84/2-' h -1 Mpc. 
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Fig. 3. — The intermittent exponent ((h) calculated by the linked-pair approximation, ( is 
the same as Fig. 2. The top scale is the comoving scale, i.e., 189.84/2- 7 h _1 Mpc. 
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Fig. 4. — The scale-dependence of Q 2n - The top scale is the comoving scale, i.e., 189. 84/2-? 
h" 1 Mpc. 
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Fig. 5. — The singular exponent 7 given by the asymptotic behavior of eq.(45). 



